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I ' Abstract 

^: 

D , We use an alternative method to the Bethe-Salpeter equation, the N-Quantum 

j^ , • approximation (NQA), for studying bound states in motion. We use this method to 

find a relativistic equation for weakly bound states of two constituents with different 
masses. We present rules for interpreting simple diagrams associated with the NQA. 
^ I We can use these rules to construct the bound state equations directly, avoiding some 

of the complications of the process. The final result is a bound state equation that 
shows Lorentz contraction in the direction of motion explicitly. This result matches 
^^ ' that of [3J found using the Bethe-Salpeter equation. We briefly discuss some other 

~ I I applications of the NQA in studying the effects of motion on bound states. 
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1 Introduction 



The purpose of this paper is to show how the N-Quantum procedure can be used to study 
bound states in motion. While our results are similar to those found using the Bethe- 
C^ ' Salpeter equation, the procedures are quite different. The effects of motion on a bound state 

have been analyzed within the framework of the Bethe-Salpeter approach in [T], [2], and 
^. While the success of the equation in describing bound state phenomena is undeniable, 
it is not without its shortcomings. As pointed out in |1], [5j, and [6j, the B-S method suffers 
from spurious solutions with negative norm amplitudes. It is also difficult to interpret the 
relative time coordinate of the B-S equation. The NQA avoids these problems and may 
also be easier to employ in certain situations. 
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The main idea of the N-Quantum approach is to expand the interacting fields that 
appear in the Lagrangian or Hamiltonian in terms of in fields. These in fields are related to 
eigenstates of the Hamiltonian with quantum numbers of freely moving asymptotic incoming 
particles. We assume that these fields form a complete set. The Haag expansion of the 
interacting fields in terms of in fields is generally an infinite series. For the purpose of 
approximation, we terminate this series, keeping only terms with a small number of in 
fields. Each term in the series contains an undetermined function of the relevant coordinates 
known as a Haag amplitude. The goal of the NQA is to derive an equation, or a set of 
equations, that can be used to solve for these amplitudes. We accomplish this by taking 
the equations of motion for the interacting fields, expanding each of the fields in normal- 
ordered products of in fields, and renormal ordering. We remove residual in fields by 
contracting with external in fields. After all in fields are contracted, the results are bound 
state equations for the amplitudes. If only low order terms are used in the Haag expansions, 
these equations will be linear in the amplitudes. 

Although they sometimes arrive at similar conclusions, the N-Quantum procedure is 
quite different than the Bethe-Salpeter method. The derivation of the Bethe-Salpeter 
equation begins with the Dyson equation for the two particle Green function, while the 
N-Quantum's roots are in Haag's operator expansion. Unlike the Bethe-Salpeter wave 
function, the position space Haag amplitudes depend only on three- vectors. With the use 
of these amplitudes, the N-Quantum procedure avoids using the relative time coordinate 
while still maintaining a covariant formalism. The Haag amplitude is similar to a Bethe- 
Salpeter amplitude with one of the constituent's mass shell singularity removed and that 
constituent's momentum restricted to the mass shell. 

The first section of this paper reviews the NQA for a relativistic bound state composed 
of two fermions in motion. A similar model for a bound state at rest composed of two 
scalars mediated by a third scalar was analyzed in [7j. Other applications of the NQA, such 
as the study of symmetry breaking, scaling limits, and a deuteron model, can be found in 
[8j, [9], and [10]. Although working in position space is in some ways more intuitive, we 
work in momentum space in this paper to make the calculations simpler. We describe the 
standard N-Quantum procedure and show how to draw and interpret diagrams associated 
with the bound state equation. These diagrams are similar to Feynman diagrams, but must 
take into account difference between off-shell and on-shell lines. We develop a set of rules 
for interpreting the diagrams to facilitate future calculations of more complicated diagrams. 

Section (3) begins with a change of variables to relative and total momentum. After 
expanding certain factors in these new variables and using some approximations proposed 
in [3], we rewrite the equation of motion to explicitly display the Lorentz contraction of 
the amplitude. Our results for the case with different masses for the two particles reduce 
to the equal mass case of |3]. We plan on discussing higher Fock state contributions in a 
future paper. 



2 The N- Quantum Procedure 

We refer to fermion 1 as an electron and fermion 2 as a proton. We begin with the equations 
of motion. The electron momentum space equation of motion (neglecting weak interactions) 
is 

(7% - mi)^i(p) = j^-r^ / dqdk5{p -k- q)Af,{q)ipi{k) (1) 

and the photon equation in Feynman gauge is 

-p^Af'ip) = —^ I dkidk25{p + ki- k2)M-ki)7^Mk2) (2) 

where we have defined ip{p) = ^p{p)^^^ and the sum over the subscript s = 1,2 is implied. 
Using Eq. (2) we can rewrite Eq. (1) as 



eV f „^,, „ S{p-k-l2 + ll)M-hh^.Mh)Mk) 



{l^Pfji - mi)V'i(p) = -7^:33 / dkdhdh (I - I )'^ ^^^ 

Since we will be substituting expansions of the interacting fields, it is important to properly 
symmetrize these fields. This is also way of imposing charge conjugation symmetry. After 
symmetrizing, the equations become 



-p'A^'ip) = ^^ J dkidk2d{p + ki - k2)[M-h),rMk2)]- (4) 

f , ^ , . ^ e'r [„„ „ S{p-k-l2 + h)[[M-h),l^Mh)]-,Mk)]+ 

(7 Pm ~ ^vVi{P) = -jTTTTs / dkdhdh 



4(27r)3 7 ^ ^ ^2-/1)2 

(5) 

We now expand the interacting fields in terms of in fields. Although only a few of the 
terms in expansions are relevant, we show some other terms that may be useful in other 
calculations as well. The Haag expansions are 



^i(p) = ^i^{p) + / d*rdH5{p + r - 6)^1 (r, b) : i)f{-r)B'''{b) : 

+ / d*rd^qd%5{p + r-q- b)J^-l'{r, b) : A^;^{q)i,f{-r)B^'^{b) : (6) 

A^'{p) = f d^rid^r25{p + n - rs) : V^r(-^i)^^(n, r2)V2"(^2) : 

+ fd%d^rid^r2S{p + ri + r2-b)^^'{ri,r2,b):tPl''{-r2)tPT{-n)B'''{b) : (7) 



where B^"'{b) is the bound state field, the colons represent normal ordering, and spinor 
indices are suppressed. The expansion for the proton field, ^2, is found by interchanging 
1 -H- 2 in the electron equation. 

Contractions between two momentum space in fields are 

< 0\rilipi)r{}}{p2)\0 > = i]/i + rmU6Mipi)d\pi +P2) 

< 0|^i;^fe)^l^(pi)|0 > = -(i/i + mi),^4-)(pi)<54(pi +P2) 
< 0\A^{pi)A'{p2)\0 > = -g^^-5i;\pi)6\pi+p2) 

where d^{p) = 6{-^p^)8{p'^ — m?). We insert the Haag expansions for the electron and proton 
fields into the right hand side of Eq. (5) and contract to find the two main terms that contain 
the factor 'ip^B^'^. (There are two additional terms that are used in renormalization.) We 
"peel" these fields off by contracting with an external ip^ and S*" field. This is equivalent 
to replacing the fields with a factor of {]/+ 1712), where p is the momentum of the proton, 
and dropping the integrals over their momenta. The result is 

RHS = Ti + T2 (8) 

where 



Ti = -y^Sm.APiMp-b + Pi) I d%,^i{l2,b)[{l^ + m2hM + ^2)? 



and 



q2^^P 

8m2ih) 
ih-pi? 



(27r)3 

^ '^mi (fc) 

(6- k-pi) 



6{p-b + pi)5rr,,{pi) / d^fc(^+mi)[(/i+m2)(7^)^2(A:,6)]^ 



2' 



Defining fi{p, h) = ^i{p, &)(]/+ m) C and using C{-j^) C = —7^, C = 7^7^, we can write 
the bound state equation in the tidier form 

flV W \f f h\ ^^ f ,i I nf^m2{p')fl{p',b) 

[J/- ]/- mi) flip, b) = / (Tp 7^^ 



2(27r)3 7 ^ ' \ fc2 

^ >lA-^+m2), (9) 

where k = p — p' and k' = b — p — p' . 

We have arrived at the desired bound state equation through a somewhat arduous 
procedure. It would be useful for future calculations to develop a process for drawing 
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Figure 1 : Graphs for the right hand side of the electron equation of motion. Heavy hues are on sheh 
and Ught hues are off sheh. The heavy double line represents the bound state (hydrogen atom). 
The empty circle represents the amplitude Fi in (a) and F2 in (b). The left fermion line is the 
electron and the right line is the proton. Similar graphs exist for the proton equation. 



and interpreting graphs, rather than contracting fields and simplifying more complicated 
expressions. The diagrams that are relevant to the preceding calculation must have one 
external off-shell electron line, an external on shell proton line, and an external on-shell 
bound state line. The two possible lowest order diagrams are shown in Figure 1. The 
rules for analyzing the diagrams associated with the N-Quantum procedure are similar to 
Feynman rules, but they must also accommodate on-shell lines. The rules are: 

1. Draw all possible relevant low order diagrams. In this case, these are diagrams with an 
"incoming" off-shell electron and on-shell proton and an "outgoing" on-shell bound 
state. Amplitude vertices should be distinct from normal vertices. Care must be 
taken to ensure the correct ordering of the following factors. 



2. Write a factor of 



for any off-shell line not connected an amplitude, where p is 



]/— m 

the momentum of the line and m is its mass (This factor is on the left hand side in 
Eq. (9)). 

3. Write ej^ for every fermion-photon vertex. 

4. Write nothing for any off shell line connected to an amplitude. 

5. Write a factor of Sm{p')fi{p', b) for the bound state vertex, where p' is the 4-momentum 
of the on shell fermion line connected to the bound state vertex and m is its mass. 
Take a transpose, then left multiply by iC and right multiply by iC if the internal 
off shell fermion line is different from the external off shell linej^ 

6. Write —gfn^/k'^ for every internal off shell photon line, where k is the momentum of 
the photon line. 



7. Integrate over the internal on shell momentum with a factor of (27r)^'^. 

8. Add a factor of (— ]/+ m) for any on shell external fermion line. 

9. Add a symmetry factor, in this case, g- 

These rules must be slightly revised for more complicated diagrams, such as those including 
off shell internal fermion lines not connected to any amplitude. 

A noticeable feature of Eq. (9) is the inclusion of both positive and negative mass shells 
on the right hand side. The opposite mass shell was dropped in |11| . a paper that uses the 
N-Quantum to study the hydrogen atom at rest. It is easy to see why the opposite mass 
shell is negligible. Looking at the fc^ in the denominator, we find when the momentum p is 
on the opposite mass shell as p' 



e 







2 2|b|g||' 



+ ^2.|1+ '^^ ^^2 ^^ |]'-(P-P^^' 



where e = -y/b^ + (mi + 7712)^, ^2 = ^"jf^ and q is the relative momentum defined in the 
following section. When the two momenta are on the same mass shell, we find 



k' 



{ki + T'kp^O{a'), 



a result found in the following section. The opposite mass shell term is suppressed by a 

2 

factor of order 0(%) relative to the other mass shell. While the negative mass shell does 
seem to be small compared to the positive shell, it may be important when analyzing higher 
order contributions. 

We can show that the second term on the RHS of Eq. (9) is equal to the first term to 
lowest order by finding a relation between the amplitudes Fi and F2. We find this relation 
by starting with the equal time anticommutator 

[Vi(x,t),V'2(y,t)]+ = o, 

Fourier transforming, Haag expanding, and contracting in fields. The relation between the 
two momentum space amplitudes is 

6mAl)Ml,b) = -5mAp)C'^fiiP,bfC. (10) 



^The full mass shell function 5m{p') is a result of symmetrization and specific to this example. For 
higher order terms, expressions will exist with higher order factors of Sin (p)- Symmetrization creates 
specific functions of the mass shell delta-functions such as Qi{pi,p2) found in Eq. (18). 



where I = b — p. Using this, the second term becomes 



' ' ^..^. JmAp')c''Mp',bfc _^^ 



2(2^)3 7 ^ ' k 

e' f,A>Jm.Ab-plfi{b-p',b) 



d P T 772 7/.(-i^+ ^2) 



2(27r)3 7 "^ ' fc' 

= ^(^ y '^ ^ ^ — p — ^^^-^^ "^^) 

This term exactly matches the first term. It should be noted that this approximation 
used only lower order terms in the Haag expansion. We could have kept the bound state 
equations as a set of coupled equations and solved them numerically for a more exact result. 

3 Hydrogen in motion 

In this section, our goal is to arrive at an equation that shows the Lorentz contraction of a 
moving bound state. This equation will be similar to that found in [3], but with different 
masses for the constituents. Although the results resemble each other, the derivations are 
quite different. We must first define the relative momentum by 

q = p-K2b (11) 

where ^i = ^ and M = mi + 771-2. In terms of the relative momentum, the proton 
has momentum q + ^26 and the electron has momentum b — p = —q + nib. We use the 

notation Ep = a /p^ + Tnf to distinguish between the energies of the two particles, and the 
subscripts || and _L to indicate components parallel and perpendicular to the bound state 
motion. 

We begin with Eq. (9) with the approximation of Eq. (10). We can left multiply by 
the inverse of the kinetic energy operator on both sides to get 



e^ (y-y+mi) 
(27r)3 (6 - p)2 - mf 



fiiP'b) = ,;_,, ;" "^J '"41 ldyrh{p',b)6mAp'h, (^) (-]/+m2). (12) 



Manipulating the denominator outside of the integral, we find 
(6 _ pf .^\ = ^E- E^f - (b - p^2 ^.2 

,2 , oA 17^, 1^(2) A o.f(2) 
r2/,, ,, \ , ^/„,3 



e^ + 2AE{e - E]^l^,^^) - 2eE]^l^,^^ _ b^ + 2b • (q + ^^b) 
+ M2(K2-Ki)+0(a3 



e^ + 2eA^(l - ks) - 2e[K2e 1 + 



_qy 2bgj| /2bg|| 



K2 y K2 



+ 2b • q - (b^ + m2)(ki - K2) + 0(a3) 






(13) 



where AE = E - e , e = Vh^TM^ , /? = ^ , and 7 = (1 - /32)"i/2. The denominator of 
the integrand is 



A;2 = fe02 _ k2 = fe02 _ k2 - kl 



We approximate /c'' as 



\2 _ 2|b|g|| 



(-) 



|b|(g|| -gy') 



2e2 



+ 0(a3 



«;2e 1 + 



/q^N2 _ 2|b|g||' 



K2 I l2 



2e2 



(/3fe|| 



Using this, we find 



fc2«-(ki+7-^fc2). 



a result quoted earlier. 

Inserting Eqs. (12) and (13) into (11) and doing the p' integral gives 



1 e^ 

(qi +7"^gn - 2KiK2AEe)fi{p,b) = 77^^(^-i/+ mi) / d'p 



K2 



(2vr)= 



(14) 



3 / 1 1 

2i^(2)ki+7-2fc2 



^l^fi{p,b)j^,{-^+m2). 



(15) 



We next follow the analysis of [3] and approximate 

(I/- j/+ 7m)rh{p', bh^{-j/+ 7712) = 442)i^(^pA+(b - p)7^/i(p', 6)7mA"(p) 

= AE('^Ei'\j-'Mp',b). (16) 

where A^(p) = — — p^'T'P "^ ^ 'Pq lowest order in a, we can write -E'^_ = ^le, E]^, = 

(2) 
Ep = K2e. Using Eq. (15), our final equation has the form 

(i;(,i +.-,|, - AM,A,P,.) = -^/.3,._W|;1|_ (,„ 

where A Af = 7A£' and ^ = ^^.^ . The Lorentz contraction in the direction of motion is 
explicit in this equation. Our result reduces to that of [3] in the equal mass case. 

4 Summary and future work 

We used the NQA to find an equation for a moving bound state consisting of two fermions 
of different masses. After some approximations, we cast this equation into a form where the 
Lorentz contraction was evident. This bound state equation matched that found through 
the Bethe-Salpeter procedure. We did not elaborate much on the final answer, i.e. explain 
how to put the Dirac structure back in or show how the frame dependence can be removed 
through rescaling variables, because such things are well discussed in |3j. The purpose of 
this paper was simply to promote an alternative method for arriving at the same answer. 
Whether the N-Quantum was the simpler process in this case is debatable, but we feel it may 
be more useful in some other calculations. Integrating over the mass shell delta- functions 
generated by the on shell lines connected to amplitudes should be easier than the alternative 
in many cases. The N-Quantum also avoids the complications of the Bethe-Salpeter method 
discussed in the introduction. 

In a subsequent paper, we will show how the N-Quantum can be used to calculate 
higher order Fock state contributions. These terms can be expressed in terms of the lower 
order amplitudes with some approximation. As an example, we can calculate the ma- 
trix element corresponding to the projection of the bound state onto the ep'j Fock state, 
M^{pi,P2,k,h) =< B\tpia{pi)tp2i3{P2)A^{k)\0 >. After using the equations of motion 
to approximate some of the higher order amplitudes in terms of lower order ones, Haag 





Figure 2: Graphs for the matrix element < -B|'0i(pi)V'2(p2)^^(fc)|O >. The left fermion line is the 
electron and the right fermion line is the proton in each diagram. 



expanding the interacting fields and contracting all in fields, the result is 



(2^)^ 



T6^ipi+P2 + k-b) 



-j/i + mi)Yfi{P2, b) [ -^ 



C'fi{pi,b)Cr{-^2 + m2) 



p\ — vri\ 
Gi(/c,pi) Gi(pi,p2) 



/fc2 



2 2 

P2-'^2 



A;2 



Ife0=£b 



C 



(18) 



where Qiih^kj) = l5mAh)Sm^{kj){l + e{k^)e{k'^) +e{-k'^)e{-k^)). The C's surround the 
transposed /2 amplitude because the electron is arbitrarily given the "a" index. The results 
could have been left in terms of higher order amplitudes for a more exact, but more difficult 
to calculate solution. We could also find this expression by interpreting the diagrams shown 
in figure (2) using the rules given in section (3). 

It would be interesting to use the N-Quantum procedure to see whether or not classical 
Lorentz contraction takes place in higher order Fock state amplitudes. This method can 
also be used to study other bound state models in which Lorentz covariance has been 
established. The subject of Lorentz covariance and bound states in motion has been studied 
in a number of papers already [13] [14] [15]. Some of the models within these works have 
interactions that lead to Lorentz invariant solutions [H], while others do not [15j. It would 



10 



be interesting to analyze some of these models in the N-quantum framework, and determine 
which interactions result in Lorentz contracting solutions. We hope that this paper has 
shown the utility of the N-quantum procedure in studying such models, and we plan to use 
it to gain a better understanding of bound state motion in future work. 
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